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Abstract. We demonstrate that the counting statistics of currents in 
periodically driven ergodic stochastic systems can show sharp changes of some 
of its properties in response to continuous changes of the driving protocol. To 
describe this effect, we introduce a new topological phase factor in evolution of the 
moment generating function which is akin to the topological geometric phase in 
evolution of a periodically driven quantum mechanical system with time-reversal 
symmetry. This phase leads to the prediction of a sign change for the difference of 
the probabilities to find even and odd number particles transferred in a stochastic 
system in response to cyclic evolution of control parameters. The driving protocols 
that lead to such a sign change should enclose specific degenerax;y points in the 
space of control parameters. The relation between topology of the paths in the 
control parameter space and the sign changes can be described in terms of the 
first Stiefel- Whitney class. 
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1. Introduction 



When a system is driven using a periodic protocol, dimensionless currents can be 
introduced that describe the fluxes per the driving period, rather than per unit time. 
The appearance of stochastic currents in response to periodic changes of parameters 
is generally referred to as the stochastic pump effect [1]. Statistical properties of 
such currents usually change continuously with changes of the shape of the parameter 
path. This is, to some degree, expected since usually small changes of kinetic rates 
in an ergodic Markovian stochastic system, which consists of only a finite number of 
discrete states, do not lead to sharp changes of the system behavior. Discrete changes 
of current characteristics can be usually achieved in stochastic systems either with 
formally infinite number of states or in the limit of low temperature [2 El when 
fluctuations are suppressed and noise cannot smear the behavior near the points of 
parameter degeneracies. In this Letter, we confront this view and demonstrate that 
some of the statistical characteristics of currents can depend discontinuously on the 
choice of the driving protocol, moreover, this property can be found in finite-size 
ergodic Markov chains. 



2. Moment generating function of currents in a 2-state model 

Consider a simple model illustrated in Fig. [l] a particle can randomly jump between 
two sites along two paths. Each path can be passed in both directions. We assume the 
system to be coupled to a heat bath at constant temperature, so that the kinetic rates, 
k" of transitions from state j to the other state via the path a, with {j,a = 1,2), 
satisfy the detailed balance condition and can be parameterized as = Kjga, where 
Kj = e^^^ , ga = e~'^^°' , with /3 = 1/(A:_bT) being the inverse temperature, whereas Ej 
and Wa are referred to as the j-th potential well energy and the a-th path potential 
barrier height, respectively. 

We will study the currents that circulate in a counterclockwise direction, see Fig[TJ 
under the action of a cyclic periodic perturbation. The currents are stochastic and, 
generally, can be characterized by two sets of probabilities TTj{n;t) that at time t the 
system is at state j, having performed n transitions during the time segment [0,t] 
through the link 1 (with the barrier Wi), counting counterclock/clockwise transition 
with the ± sign, respectively. Hereafter, we will notationally suppress the time 
parameter t in Trj{n;t) and other variables. The Master Equation for TTj{n) is given 

by 

7ri(n) = -{kl + kf)ni{n) + fc^7r2(n - 1) + k^n2{n), 

(1) 

Mn) = -{H + fci)7r2(n) + kliriin + 1) + fc?7ri(n). 
Upon a Fourier transform with respect to n, Eq. ([T]) translates into 

oo 

Z = H{x;t)Z, Z = {Z,,Z,), Z,{x;t)= ^ 7r,(n; i)e'^", (2) 

n— — CK 



with the "twisted" evolution operator 

fy^^^ +ff2) K2(5ie'^ +.92) A .o^ 

~ \ '^ligie-''' + 92) ~n2{gi+g2) J' ^' 

The averaged over the final state MGF Z{x; t) = Zi{x; t) + ^2(x; ^) that contains 
complete information on the current distribution function, including the statistics of 




W2 

Figure 1. A model of stochastic transitions between two states along two possible 
paths. Different paths are characterized by the values of the corresponding 
potential barriers Wi and W2- The two states are characterized by the values of 
their well depths, Ei and £2- Periodic modulation of these parameters leads to 
appearance of a nonzero, on average, circulating current in the preferred (clockwise 
or counterclockwise) direction. 

rare events, is obtained by solving Eq. ([2| and can be expressed in terms of a time- 
ordered exponential [5] 

Zix;t) = (l|Texp (^^ Hix,t')dt'^ k), (4) 

where (1| = (1, 1), and |7r) = (tti, 1:2) is the vector of the initial system populations. 

3. Degeneracy point of 7?(x)-eigenvalues 

The lowest current cumulants are described by the properties of Z{x) near the x — 
point, and their behavior in the limit of adiabatically slow evolution of parameters 
is always defined by the eigenvalue of H{x) with the largest real part and the 
corresponding eigenvector [B]. In the vicinity of x = 0, the eigenvalues of H{x) 
are nondegenerate, provided all kinetic rates in ([s]) are nonzero (finite temperature). 
This follows from an observation that H{0) describes the evolution of probabilities in 
a 2-state Markov chain [6] . Such an evolution matrix has a unique zero-mode solution 
corresponding to the unique steady state of an ergodic Markov chain. This uniqueness 
guarantees that the eigenvalue of H{x) with the largest real part is nondegenerate for 
sufficiently small x- 

We will study the effect of geometric phases on the full counting statistics of 
currents that does not apparently appear when only lowest cumulants of the current 
distribution are measured, namely, the matrix ([S]) has a unique point in the parameter 
space, where its eigenvalues are degenerate. This happens when simultaneously three 
conditions x = tt, Wi = W2, and Ei = E2 are satisfied, which can be seen by explicitly 
inspecting the eigenvalue problem for the 2x2 matrix given by Eq. ([3|. 

Interestingly, the degeneracy conditions do not depend on the inverse temperature 
p. The fact that the degeneracy point is encountered at a finite value of the counting 
parameter, x — '^^ means that it influences the properties of the full counting statistics, 
rather than the lowest cumulants of the current distribution. We note also that at 
X = TT the matrix H^n) has all real entries. 
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The MGF at x = tt, 

2; OO 2; oo 

Z{7:)=Po-P,= J2 ^ji^k)- J2 + (5) 

J — 1; fc— — oo J — 1; k— — oo 

has a simple physical meaning, namely, it is the difference of the probabilities Pq and 
Pi to make an even and odd number of transitions through a given link, respectively. 

The time evolution of Pq—Pi must be influenced by the presence of the degeneracy 
point. The evolution operator H{t:) has two nonpositive eigenvalues 

£± = -K+g+ ± ^K\g^_ -5-), (6) 

with K± — (ki ± K2)/2, g± — gi ± g2. Near their degeneracy point, we parameterize 
g^n+ — rsin0 and K-g+ — rcos0 which yields e± ~ — ± r for the 

eigenvalues in the limit g- ^ g+, k_ <C k+. The eigenvalues as functions of 
(r, 0) lie on a two-cone surface. Taking the higher eigenvalue £-|-, the eigenstate 
|M+(7r)) = (ui,U2), determined by solving a 2 x 2 linear problem, has a form 
{ui,U2) = a((/))(sin((/)/2), — cos((/)/2)) with a((/)) being an arbitrary periodic function 
of (j). The requirement on the eigenstate not to be zero imposes that a{4>) is either 
strictly positive or strictly negative function, however, the functions sm{(j)/2) and 
cos(0/2) change the sign when changing </> around a full cycle. If one chooses the 
smooth field of eigenstates of H{tt) over the space of all possible values of parameters 
Ei and Wj, except the degeneracy point, such a vector field has to be double- valued 
and, along a nonintersecting path around the degeneracy point, the eigenvectors of 
H(Tr) change the sign. To understand, how this effect can influence the stochastic 
currents, we should look at the evolution of MGF in response to externally induced 
periodic changes of kinetic rates in the model. 

4. Evolution of MGF 

4.I. The case of constant parameters 

When all parameters in the evolution operator are time-independent the time- 
ordered exponential in Q becomes just a matrix exponent, which can be simplified 
by introducing the left and right eigenstates {u±\ and \u±) of the matrix Ht^, 
corresponding to the eigenvalues in Then Q gives us 

Po - Pi = C+e^+* + C_e^-*, C± = {l\u±{t = 0)){u±{t = 0)|7r). (7) 

Since both eigenvalues e± are negative, the absolute value of Pq — Pi decays 
exponentially with time. Moreover, since e+ > £-, the term C-e"^"* quickly becomes 
exponentially suppressed in compared to C+e^+*. The coefficients C± are totally 
determined by the initial probability distribution and kinetic rates. Their values 
describe the initial relaxation processes to the steady state distribution and remain 
the same even if kinetic rates are changing with time adiabatically [71 . At time scales 
larger than !/(£+ — £-) the sign of Pq — Pi is determined by the sign of C+, which 
can be found and analyzed explicitly. Depending on the initial choice of parameters, 
including initial state probabilities, can be either negative or positive, however, 
the sign of Pq — Pi is fixed after the term with coefficient C_ becomes exponentially 
suppressed. In other words, the sign of Pq — Pi is the invariant of the evolution which 
can take values in a discrete two-state set. 
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Figure 2. Evolution of the parity probability difference in a driven 2-state 
stochastic system. The choice of parameters is k+ = 2, <;+ = 1, k— = 
0.8 + 2cos(27rt/T)/3, g_ = 0.4 + sin(27rt/T)/3, T = 30. The contour in the 
space of control parameters does not enclose the degeneracy point. 



4-2. Geometric phase 

We further ask a question whether the sign of Pq — Pi can be switched by adiabatically 
slow evolution of parameters g- and K- around a cycle. In the adiabatic limit, a 
naive strict quasi-steady-state approximation predicts that a slow perturbation merely 
modulates the decay exponents in Q and the vector Z{tt) is given by instantaneous 
eigenvector This means that, in principle, the sign of Pq — Pi can change during 

the evolution since the components of change with time and their sum can switch 
the sign. However, if the parameters complete a cycle, the initial and final eigenstates 
coincide, while the coefhcient C+ is time-independent. Hence, one may not expect to 
find a sign change of Pq — Pi at the end of the cycle. This analysis, however, does not 
include the effect of geometric phases, which can be non-perturbative [8|. 

First, consider the case when a cyclic path in the space of control parameters does 
not enclose the point g^ = — 0. We can always choose a smooth parametrization 
of the eigenvectors of H{tt) in the region Sc inside the contour. According to 
[H El ini Uni HI]; -fo — Pi, after periodic and adiabatically slow evolution of the 
parameters around the contour C during time T, is given by 



Po-Pi= C+ exp y J e+{t)dt + (^goom(C) j , (8) 

where the geometric contribution to the exponent in ([s]) is given by 

</5gcom(C) = / dg-dK-Fg^^^_, (9) 
Jsr- 



'Sc 

where = {dg_u+\di^_u+) - {df,_u+\dg_u+) or explicitly 

K-.9-^+.g4 

2{Klgl+n^_{gl-g^_)Y/^- 

The eigenvalues e+{t) and the geometric phase ^PgcamiC) in ([s]) are both real, 
and therefore the exponential term does not affect the sign of Pq ^ Pi- Fig. [2] shows 
the results of our numerical solution of the evolution equation ([2| at x = tt. The fact 
that Pq — Pi has the same sign during the whole evolution is not generic. Generally, it 
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Figure 3. Evolution of the probability difference Pg — Pi in a driven 2-state 
stochastic system when a contour in the space of control parameters encloses 
the degeneracy point. The choice of parameters is k+ = !,(;+ = 1, k_ = 
2cos(27rt/T)/3, g- = sm{2nt/T) /3, T = 30. 

can change but the fact that it does not change after each period of driving protocol 
is always confirmed if the contour does not enclose the degeneracy point. In other 
words, the sign of Pq — Pi is an invariant of an adiabatic change of parameters along 
a closed contour that does not enclose the degeneracy point. 

4-3. Robust switching of the sign of Pq — Pi 

We now consider the case when the contour is cyclic and encloses the degeneracy 
point (7_ = K_ =0. The previous arguments cannot be applied directly since at the 
degeneracy point the adiabatic approximation breaks down and the parametrization 
of the eigenstates inside the contour cannot be chosen smoothly. They do, however, 
imply that sgn(Po — Pi) does not change upon deformations of C that do not encounter 
the degeneracy point reducing our analysis to the contours, restricted to a small 
neighborhood of the degeneracy point. Our analysis, in section 3, of the vicinity 
of the crossing point showed that the smooth choice of eigenstates along a contour 
requires that, upon completion of the cycle, the eigenvector changes sign. Physically, 
this change of sign means that the probability difference Pq — Pi changes the sign. 
Our numerical solution of Eqs. ^ for evolution along a contour that encloses the 
degeneracy point, is shown in Fig. |3] It confirms the sign change prediction. Note 
that this result is nonperturbative. It cannot appear as an effect of small non- 
adiabatic corrections to the decaying exponent in ([8|. The factor —1 is a topological 
factor rather than a geometric phase effect in a sense of Ref. [6j. One can see 
from Fig. |3] that this factor does appear after an adiabatic cyclic evolution of the 
control parameters that enclose the degeneracy point. More generally this implies that 
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sgn{Po—Pi) — (— l)"('^^sgnC+ with n{C) being the number of times the cycle C winds 
around the degeneracy point. This property provides an approach to manipulate this 
sign with high robustness against perturbations of a driving protocol. 

5. Topological interpretation of the invariant and generalizations 

The sign of the probability difference in ([5| is, in fact, a topological invariant of the 
driving protocol, the latter being viewed as a cycle in the space M = \ {0} of 
parameters x = {g-,K-) where degeneracy of H{'K]x) never occurs. By associating 
with X & M the ground state of i?(7r; a;), the former defined up to a real multiplicative 
factor, we build a linear real fiber bundle E over M, or equivalently a map g : M ^ 
MF^. Here the projective space EP^ is interpreted as a circle of "normalized states" 
+ (^2)^ = 1 with the opposite points on the circle identified, since there is no 
preferred way to define the sign of the ground state. The map g associates with a 
driving protocol C a cycle g{C) in MP^ = , and each time C winds around the 
degeneracy point, g{C) produces a single rotation in MP^, so that n{C) = n{g{C)). 
The last statement can be reformulated in terms of the first Stiefel- Whitney class 
Ci = Ci{E) [12] of E that can be viewed as a Z2-valued function on cycles in M, so 
that Ci(C) = n{g{C)) mod 2 and the invariant {-l)"-^^) = (-l)'=i(C) of the counting 
statistics is obtained by evaluating the Stiefel- Whitney class at the driving protocol. 

In this letter we studied in some detail a Z2-invariant in a simple 2-node system. 
In a general network represented by a finite connected graph there is a set of invariants, 
we are going to describe briefly here, with the analysis of the general case postponed 
to a further publication. With any "independent" set a. = {ai, . . . ,ak] of links (i.e., 
if they are all withdrawn the network still remains connected) we can associate a Z2- 
valued observable /^(C) = sgn^^(— °jPo(a; C), where a = {ai, . . . , a^} is a set 
of binary variables, and Pa{cy.; C) is the joint probability that for protocol C the total 
number of jumps over link aj was even/odd for Oj = 0, 1 respectively. The protocol 
is assumed to go through the space Mq, of parameters that avoid the ground-state 
degeneracy of the twisted operator with X7 = for chosen links 760;, and 

X'y = 0, otherwise. Similar to the simple two-node case the ground state defines a linear 
bundle Ea over Ma with the first Stiefel- Whitney class ci — ci{Ea), or equivalently 
a map g : Mq, — > MP^"^ with N being the number of nodes in the network. Viewing 
RP^~^ as an {N — l)-sphere with opposite points identified we have a two- fold cover 
5*^^^ MP^^^. Since all cycles in S^~^ are contractible, there are two equivalence 
classes of cycles s in the projective space: the ones produced by cycles in S'^~^ and 
the ones produced by contours that connect opposite points in S^^^; this is referred 
to as [s] = and [s] = 1, respectively. The topological factor associated with a driving 
protocol is given by (— l)!^^^)] = (_i)ci(C) ^hat the counting-statistics invariant 
Ia{C) = (— l)'^'^^*-^^sgnC+, with C+ depending on the initial distribution only, being 
independent of the driving protocol. 

The topological factor is very similar to the topological Berry phase in quantum 
mechanics of systems with time-reversal symmetry (no magnetic fields), where H{x) 
is a family of Hamiltonians, represented by real Hermitian operators [T51 15]. In our 
case the operators are real, however, apparently non-Hermitian. They are in fact 
Hermitian with respect to a scalar product that depends on the kinetic rates. The 
latter results in the geometric contribution </3gcom(C') in Eq. (|8| that does not affect 
sgn(Po-Pi). 
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6. Discussion 

We demonstrated that some of the properties of the full counting statistics of currents 
in finite size stochastic models can be controlled by periodic changes of the control 
parameters without continuous dependence of the result of operation on the choice of 
the path in the parameter space. This finding reveals the importance of the degeneracy 
points of the eigenvalues of H{x)j governing the evolution of the current MGF. 
Although the effect is the property of the full counting statistics [3 [HI [151 [El [13 [E] , 
it is measurable experimentally, since the probabilities of rare events in the considered 
model are suppressed at least exponentially. As a result, only a finite number of 
terms in the definition of the parity difference ([s]) should be determined, which is 
achievable by repeating the measurement of the number of transition through a link 
after sufRciently many cycles. The simplicity of our model suggests that the topological 
geometric factors could be a generic property of current statistics. The phenomena 
that they represent in general situations of multistate or Langevin evolution can be 
much more complex than the sign switching that we demonstrated. 
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